How to De-Rotate the Cosmic Microwave Background Polarization 
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If the linear polarization of the cosmic microwave background (CMB) is rotated in a frequency- 
independent manner as it propagates from the surface of last scatter, it may introduce a B-mode 
polarization. Here I show that measurement of higher-order TE, EE, EB, and TB correlations 
induced by this rotation can be used to reconstruct the rotation angle as a function of position 
on the sky. This technique can be used to distinguish primordial B modes from those induced 
by rotation. The effects of rotation can be distinguished geometrically from similar effects due to 
cosmic shear. 
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The polarization of the cosmic microwave background 
(CMB) can be decomposed into gradient (E mode) and 
curl (B mode) components [lj. Primordial density per- 
turbations produce a polarization pattern that is purely 
E mode at the surface of last scatter, while primordial 
gravitational waves, such as those from inflation, produce 
a B mode [2j . There is now an active experimental effort 
to detect B modes, as these would constitute a "smoking 
gun" for inflationary gravitational waves Q . 

There may, however, be other mechanisms — apart 
from gravitational waves — for producing B modes. The 
most widely considered is cosmic shear: the deflection 
of CMB photons due to weak gravitational lensing by 
density perturbations along the line of sight will convert 
some of the E mode at the surface of last scatter to a 
B mode [H . Cosmic shear of the CMB is no longer the 
realm of futuristic theorists; it has recently been detected 

Another possibility is a rotation of the linear polar- 
ization of the CMB as it travels from the surface of last 
scatter [5]. This could occur, for example, if there is 
a quintessence field that couples to the pseudoscalar of 
electromagnetism [6j. In this case, the polarization would 
be rotated by an angle a that is uniform across the sky. 
A fraction sin 2a of the E mode would thus be converted 
into a B mode [5[. This B mode could be distinguished 
from that due to gravitational waves by the parity break- 
ing EB and TB cross-correlations that it produces, as 
well as by the shape of the TB, EB, and BB power spec- 
tra. Searches with current data already constrain such a 
uniform rotation to be no more than a few degrees 0, S] • 

But what if the rotation angle a varies across the sky? 
Refs. fioj have recently proposed models in which this 
might occur by virtue of a spatially and time-varying 
scalar field coupled to the pseudoscalar of electromag- 
netism, and something similar may arise from dark- 
matter magnetic moments (Faraday rotation could 
also rotate the polarization [9( , but this can be identified 
with multifrcquency maps.) If the distribution of rota- 
tion angles is symmetric about zero, then there will be a 
B mode induced with no parity-breaking TB or EB cor- 



relations [lOj . Can this type of B mode be distinguished 
from that due to primordial gravitational waves? 

In this Letter, I show that a position-dependent ro- 
tation of the polarization induces higher-order correla- 
tions in the temperature-polarization of the CMB. I then 
show how these new correlations can be used to measure 
the rotation angle a(n) as a function of position n on 
the sky. The observed polarization pattern can then be 
corrected to construct the primordial polarization pat- 
tern. It will thus be possible to distinguish whether a 
B-mode pattern, if detected, is primordial or due to a 
post-recombination rotation. The techniques I will dis- 
cuss can also be used to test an experiment for systematic 
artifacts. 

The reconstruction algorithm will be similar to that 
developed to reconstruct the cosmic-shear field and 
also to that developed to test for statistical isotropy (SI) 
[l3l | . In fact, the effects of rotation may be viewed as a 
possible contaminant for cosmic-shear maps. However, 
the effects of cosmic shear (or SI violations) and rotation 
can be distinguished, as (a) rotation has a different ge- 
ometric effect on the polarization, and (b) cosmic shear 
affects the temperature as well. 

The linear polarization at each point h on the sky is 
quantified by Stokes parameters Q(n) and U(h), mea- 
sured with respect to the 8-cf) axes. These Stokes pa- 
rameters are components of a symmetric trace-free 2x2 
polarization tensor field V a bin) which can be expanded 
in terms of tensor spherical harmonics as, 



P*6(ft) = E E K^(L )a6 («) + B lm Y* m)ab (n) 
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Here ' F (L)a&(") and Y (lrn)a b y 
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(h) are complete sets of basis 
functions for the gradient (E mode) and curl (B mode) 
components of the polarization. Orthonormality of the 
basis functions allows us to write the expansion coeffi- 
cients as 



Elm — 



dnV ab (fi)Y^*(n), 
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Bim = J dnV ab (n)Y^*(h). 



(3) 



Suppose now that the polarization pattern at the sur- 
face of last scatter has no B mode and that the polariza- 
tion at each point h is rotated by an angle a(n). In that 
case, the observed Stokes parameters will be 
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The concordance of measured TT, TE, and EE power 
spectra suggest that a is relatively small, and that the 
uniform component of a is small [7[. We therefore 
work to linear order in a. The change to the polariza- 
tion tensor induced by the rotation is then 5V a b(n) = 
2a(n)'P^ b (n), where V T ab (n) is rotated from the original 
polarization V a b by 45°. Thus, if V a b is pure E mode, 
then V x ah is pure B mode and can be written as 



;)ab' 
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To calculate the curl component induced by rotation, 
we insert 5V a b into the expression, Eq. |3]), for Bi m . We 
then expand the rotation angle in terms of (scalar) spher- 
ical harmonics Y {lm) as a(n) = J2lm a LMY( LM) (n), to 
obtain 

Bim = ^2^2 a LM E hm2 J dnY^ * Y (LM) Y^ 3m ^ ab 

LM 12^712 

= 2(~l) m y^ QjLM^ 2 w 2 ^tof 2 ro2- g H 2 i ( 6 ) 

where the sum is taken only over li values that satisfy 
L + I + l 2 =even, 



L = ( I V L\( I V L 
"' "^-2 2 /1^0 

in terms of Wigner-3j symbols, and 



(7) 



£&# ro ,= / dnY { * lm) (n)Y {l , 



(n)r (LM) (n). (8) 



We thus see that rotation induces a B mode [H, [l(| • If 
there is a power spectrum for aLM, then the CMB power 
spectrum C BB can be calculated from Eq. ^ (loj . 

There is also an 0(a) change in the E mode induced by 
rotation with precisely the same form as Eq. (J6]) , but with 
contributions only from L + l + l2=odd, rather than even. 
In the discussion below, we focus for brevity and clarity 
on the induced EB and TB correlations. However, the 
entire discussion applies (with differences to be pointed 
out below) to induced EE and TE correlations. 



The next step is to consider the correlation of the in- 
duced B mode with the original E mode, as well as with 
the temperature T. Recalling that the original E modes 
have (assuming Gaussian initial conditions) expectation 



values (Et m E Vm , : 



= C^ E S w d mm ' (where Cf E is the EE 
power spectrum), the correlation between the observed 
E and B modes is 

(Bi m E*, m ,) = 2—^=Cf F '5u'5 mm ' 
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I have split off the L — term in Eq. ([9]) to show that 
the result for a uniform rotation angle a is (Bi m El, m ,) = 
2aCf E 8 U '8 mm i, as it should be @. The sum in Eq. © 
is taken only over L + I + Z'=evcn. The expression for 
the induced correlations (Ei m El, m ,) is the same as that 
in Eq. ©, but with L + I + Z'=odd, and no L = con- 
tribution. 

Likewise, given that the temperature T and E-modc 
polarization are correlated at the surface of last scatter, 
there will also be a nonzero TB correlation induced by 
rotation. The expression for (Ti m B*, m ,) is identical to 
Eq. ^ with the replacement C EE — > C™, again for 
L + 1 + /'=even. There will also be TE correlations of the 
same form, but with L + I + /'=odd and no L — term. 

If a(h) varies with n, then a^M ^ for L > 1, and if 
so, there will be correlations between Ei m and By m i (and 
Tim and Bv m ') of different Im and I'm' . The existence 
of these off-diagonal correlations can be used to measure 
each of the rotation multipole moments a/ m , and thus 
a(n). The relevant formalism is similar to that for mea- 
suring the cosmic-shear field or for searching for SI 
violations 13(, so we can adopt results from prior work. 
To do so, we note that Eq. §§§ is identical to Eq. (Al) 
in Ref. [l3( with the identification X = B (or T in place 



EE 



(orC, 



TE 



of B), X' = E, and D^,' = a LM C\ 
place of C EE for TB). 

Our goal is to obtain the minimum- variance estimator 
aLM that can be obtained from a full-sky polarization 
map or a temperature-polarization map. We suppose 
that the maps are provided as a measured temperature 
T map (n) and Stokes parameters Q map (n) and C/ map (n) 
in A^pix pixels on the sky. The temperature (polariza- 
tion) in each pixel receives contributions from the sig- 
nal, which is the temperature (polarization) on the sky 
smoothed with a Gaussian beam of full-width half max- 
imum (FWHM) 0f w hmi and a Gaussian noise with vari- 
ance dip (<Tp). The power spectra for the map are then 

q ,map _ | ]/(/ 2 | £i A _|_ £^ w J iere (J > n j g fae no j se power 

spectrum for A (e.g., A = {TT, EE, BB, TE, TB, EB}). 
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TB,n 
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EB.n 
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smearing is taken into account with the window func- 



tion Wi = exp(— l 2 a 2 /2), with cr b = O iwhm /VS In 2 = 
0.0742 (0 fwhm /l°). 

We now derive the minimum-variance estimator oilm 
that can be obtained from EB correlations; the results for 
the estimator that can be obtained from TB correlations 
will be identical with the replacement E^B; TE and EE 
estimators are similarly derived. 

Following Ref. (l3| . the minimum- variance estimator 
for each D W M = 2a L MCf E H w that can be obtained 
from the polarization map is (see also Refs. [l5| ) 



Ar M ' map = (G^r 1 E B? m ap (E™vy& m „ (10) 



2Cf E H W . The second is 



where Fh 



LM 



(14) 

where the matrix inversion is the in the 2x2 EB-BE 
basis. Note that the superscripts A on D are necessary 
in these expressions, as these quantities differ for EB and 
BE. The variance to the first estimator is 
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with 

4tt(2L + l)Gf v = (21 + l)(2l' + 1) (C^of , (11) 

where Cji££, m , is a Clebsch-Gordan coefficient. 

The coefficient (Xlm can be estimated from mea- 
surement of -Dy/ M from each IV pair through a^M — 
D W M I (2Cf E i/ z f,). One can then average the estimates 
of aw from all of the IV pairs. The trick, though, is to 
weight these all in a manner that minimizes the variance 
to oilm- If each estimator £)^ M ' ma P we re statistically in- 
dependent, then we could simply weight by the inverse 
variance. However, things are a bit (though not much) 
more complicated. 

Each of the estimators £)^ ( M ' ma P are statistically in- 
dependent for different LM. They are also statistically 
independent for different IV , except that £)^ ; M > ma P [ s cor _ 

related with Df J ,f I ' map . To take this into account, we take 

V > I and then consider for I ^ V EB modes as well as 
BE modes. For I — V, there is only a single variance; for 

V > I, there is a 2 x 2 covariance matrix in the EB-BE 
space. 

Write the covariances between the different D w ' map 



nil' 
L AA' 



Ggi / T;iM,A,mapMM,A',map\ r-,r,\ 



for {A, A'} = {EB, BE}. For I = V, there is no dis- 
tinction between EB and BE; the variance in this case 



is then C^ EjBE 
covariances are C 
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For V > I, the 
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We now write two estimators, a 1 ^ and a^ f , the first 
coming from EB correlations with I = V and the second 
from those with I < V . We will then average them, with 
inverse-variance weighting, to obtain the final estimator. 
The first is 



?LffXi \2 p ) LM,ma,p r ,L I nil 

( - r ii/ L 'BE,BE 



u 

BE, BE 



(13) 



The variance (cr i<i> ) 2 to the second estimator is obtained 

01 LM 

from 



AA' 



AA' 

(16) 

The final minimum- variance estimator ulm is then ob- 
tained by averaging, with inverse-variance weighting, the 
two estimators above: 



-i=V 

LM 



OILM 



(^)" 2 + s S(V<'') 



i<r a '=i') 2 + (v a ><i'Y 



(17) 



The variance (cr aLM ) 2 for this estimator is given by 
k ul ) _! = (v a '=''y 2 + (v«')~ 2 - For small L (e.g., 

LM LM 

a rotation dipole), the two terms will contribute compa- 
rably to the statistical weight. For larger L, the V > I 
estimator should carry most of the statistical weight. 

Estimators from TB correlations are identical to the 
EB estimators discussed above with the replacement 
E^B. Likewise, there will be EE and TE correlations 
similarly induced. (Things simplify for EE, as the 2x2 
AA' covariance matrix becomes a single variance.) The 
estimators for aLM for EE and TE and their variances 
can be constructed analogously. There will also be BB 
correlations, but they will be higher order in a. There 
will be no TT correlations induced, as the rotation does 
not act on the temperature. Since T and E are corre- 
lated in the primordial polarization field, there will be 
cross correlations between the estimators o~lm from EB, 
TB, EE, and TE. This may be an order-unity effect if the 
statistical weights of the various estimators are compa- 
rable. The expressions for the complete covariances are 
long, and so I leave them for future work. 

Here I have shown that CMB temperature-polarization 
statistics can be developed to measure the angle a(n) 
by which CMB photons were rotated, en route from the 
surface of last scatter, as a function of position n on the 
sky. Explicit formulas to obtain the coefficients ollm i n 
a spherical-harmonic expansion of a(n) from a full-sky 
CMB map were provided. This technique can then deter- 
mine whether B modes, if detected, occur at the surface 
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of last scatter or are due to a post-recombination rota- 
tion of the polarization. It is interesting to know that the 
rotation angle can be determined from the data, rather 
than by assumption. And if the rotation angle is assumed 
to be zero, then the techniques developed here can pro- 
vide a test for systematic artifacts in the data. (Ref. (lt3 | 
suggested tests for systematics along these lines.) 

The rotation and cosmic-shear formalisms share some 
similarities, and so rotation, if it exists, could show up 
as an artifact in a cosmic-shear analysis. However, the 
detailed effects are different and can be distinguished in 
the data. First of all, cosmic shear has a different parity 
than rotation; a given L mode of the cosmic-shear field 
correlates I and I' modes of E and B, respectively, only 
if I + I' + L is odd, while rotation correlates them only 
for I + V + L odd. Furthermore, cosmic shear acts on 
temperature and polarization, while rotation leaves the 
temperature map unaltered. 

The formulas for the estimators Slm will need to be 
modified to take into account partial-sky coverage in a re- 
alistic map. However, it will be straightforward to adapt 
the techniques that have been developed to measure cos- 
mic shear of the CMB on a partial sky to measure the 
rotation angle. Likewise, it is straightforward to simplify 
the full-sky analysis performed here to the flat-sky limit, 
which may be appropriate for sub-orbital experiments 
that map the CMB on a small patch of sky. Since the for- 
malism to reconstruct the rotation angle resembles that 
to determine the cosmic-shear field, there may be other 
cosmic-shear techniques that can be adapted for rota- 
tion. For example, maximum-likelihood techniques [HI 
may be developed to provide even more sensitive probes 
of rotation than the quadratic estimators discussed here. 

Finally, it is of interest to know quantitatively how 
well the estimators presented here can be used to con- 
struct the rotation angle. Evaluating the expressions for 
<Ja LM for WMAP values for up and 9f w hm, we find for 
the TB estimator u aLM = 7° for L = 2 [19j, consistent 
with the current WMAP la constraint to a uniform rota- 
tion (recalling that «la/ = V^na for L = 0). This then 
increases, by about 50%, to L — 100. At higher L, the 
noise increases due to WMAP's finite angular resolution. 
For WMAP, the values for a aLM for the EB estimator are 
seven times larger, and thus not constraining. Using val- 
ues for up and 9i w hm appropriate for the Planck satellite, 
we find that the errors to the EB and TB estimators are 
more comparable; e.g., a aLM = 1.6° and 1.3° for EB and 
TB, respectively, for L = 2 rising slowly to cr aLM = 3.4° 
and 1.9° for L = 500. The noise then increases rapidly 
for / > 500, when the correlation angle becomes smaller 
than the polarization correlation angle. More detailed 
and comprehensive numerical results will be presented in 
Ref. 0. 
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